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Abstract. Spectroscopic studies of Planetary Nebulae (PNe) and H II regions have driven 
much development in atomic physics. In the last few years the combination of a generation of 
powerful observatories, the development of ever more sophisticated spectral modeling codes, and 
large efforts on mass production of high quality atomic data have led to important progress in 
our understanding of the atomic spectra of such astronomical objects. In this paper I review such 
progress, including evaluations of atomic data by comparisons with nebular spectra, detection 
of spectral lines from most iron-peak elements and n-capture elements, observations of hyperfine 
emission lines and analysis of isotopic abundances, fluorescent processes, and new techniques for 
diagnosing physical conditions based on recombination spectra. The review is directed toward 
atomic physicists and spectroscopists trying to establish the current status of the atomic data 
and models and to know the main standing issues. 



1. Introduction 

Nearly all modern observational research in gaseous nebulae (PNe, H n regions, and circumstellar 
nebulae) involves some kind of spectra, whose interpretation requires some understanding of the 
atomic processes involved and the handling of atomic data. The level of understanding of 
such processes and the data required for various applications ranges from phenomenological 
and statistical studies of commonly prominent lines, to compilation of raw atomic data for 
direct inspection (such as line finding lists), and to modeling of synthetic spectra attempting 
to fit the observed spectra. On the other hand, the accuracy and quantity of atomic models 
describing various features in observed spectra have evolved with time, driven by great advances 
in ground- and space-based observatories, computer technology and experimental techniques. At 
present there are models and data accurate enough for satisfactory analysis of at least the most 
prominent spectral features. Also the advent of on-line databases and spectroscopic tools has 
revolutionized the dissemination of results of atomic physics research. However, the demands 
on the quantity and quality of the atomic data will grow as new instruments delivering greater 
sensitivity and spectral resolution become available. 

A review of atomic processes spans an audience that includes specialists in the study of 
atomic systems with tools capable of providing models and atomic data, nebular astrophysicists, 
astronomers seeking to understand the reliability and accuracy of their modeled spectra, and 
those who use the models and data to compare synthetic spectra with observations, diagnose 
the physical conditions and compute chemical abundances of nebulae. The present review is 
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aimed at researchers in atomic physics and spectroscopy who wish to understand the reach of 
currently available atomic models and identify present needs. 

The physical conditions in gaseous nebulae of interest here are those of photoionized plasmas 
as observed in the ultraviolet (UV), optical, and infrared (IR) bands. Roughly speaking these 
conditions are: temperatures of the order of 10 4 K and electron densities between 10 3 and 
10 8 cm -3 . This review will not treat the X-ray band which in gaseous nebulae results from 
mechanically heated coronal gas. A good review on this subject can be found in pQ. 

The present review is organized according the nuclear charge of atomic species, starting with 
hydrogen and helium, and following with the second and third row elements, the iron-peak ions, 
and then onto the heavier elements. I then finish the paper with a discussion of various general 
issues of interest. 

2. Hydrogen 

The spectrum of hydrogen in H n regions has been studied extensively for quite some time. 
In principle the spectrum is determined by the recombination rates into each level and the 
subsequent radiative cascades. Such a process is accurately modeled in two idealistic situations, 
the so called "case A" in which all lines are optically thin and "case B" which assumes that the 
optical depth of Lya goes to infinity. Two complications arise beyond these approximations: 
(1) a detailed solution to the radiative transfer of Lya photons including the effects of self- 
absorption, removal and collisionally induced transitions between the 2s and 2p states, and (2) 
when the electron temperature and density of the plasma are high enough for effective collisional 
excitations on to n > 2 states. 

Dennison et al. [2] discuss the first of these problems and propose an observational test. They 
claim that radio-observatories will soon be able to detect the 2si/ 2 — %Pi/2 and 2s\/2 — 2p^/ 2 
lines at 1.1 GHz and 9.9 GHz respectively. Dennison et al. explain that removal of Lya photons 
by dust would limit the pumping of H atoms into the 2p states, and under these conditions the 
2si/2 ~~ transition will appear in stimulated emission while the 2si/ 2 — ^Ps/2 nne wm appear 
in absorption. In general, the relative strengths of these two lines should serve as diagnostics of 
the populations of the 2s and 2p states and the responsible processes. Further, Dennison et al. 
suggest that removal by dust is the dominant mechanism acting on Lya photons in H II regions, 
which if correct would circumvent the need to solve the radiative transfer problem. 

In regards to collisional excitation of hydrogen, Pequignot & Tsamis [3] study the evolution 
of calculated collision strengths for Is — > n = 3,4, 5. They point out large variations in the 
theoretical values up to the last calculation of Anderson et al. [4 J . Pequignot & Tsamis compare 
their observations for the PN G135.9+55.9 with predictions of models using various collisional 
data sets. They find that the collision strengths of Anderson et al. yield the most consistent 
results, yet these may have not reached ultimate accuracy. 

Modern high-signal-to-noise optical spectra allow for accurate measurements of lines and 
continuum around the Balmer and Paschen recombination series. By fitting the results of 
detailed spectral models to these measurements it is possible to determine the electron density 
and temperature in the nebula [5J. Surprisingly, temperatures obtained from the Balmer series 
often disagree with the temperatures derived from line ratios of collisionally excited lines, such 
as the ratio of nebular to auroral lines of [O in] . Moreover, for a large sample of PNe studied by 
Zhang et al., the hydrogen Balmer temperatures are typically lower than the collisional oxygen 
temperature by several thousand degrees. These temperature differences are used in support 
of the so called 'temperature fluctuations' or 'temperature variations' in nebulae. However, a 
flag of warning on the use of hydrogen temperatures comes from the fact that in four PNe 
where temperatures from the Paschen series were determined, these disagree with the Balmer 
temperatures at the three sigma level or more. 



3. Helium 

Modeling the He I recombination spectrum has become increasingly reliable, with an accuracy on 
line emissivities of the order of a few percent (see [SI [8]). Yet, this is not at the 1% accuracy 
level needed to place useful constrains on the primordial helium abundance in cosmological 
studies [HI [10]. Recently, Bauman et al. [6] worked out the recombination problem in fine 
structure in the low density limit. They found no significant effects on the recombination 
spectrum from spin-orbit coupling. However, they found that some of lines may be uncertain 
due to use of poor quality atomic data, specifically the photoionization cross sections for states 
with 9 < n < 20 and L < 3. 

A practical test on the accuracy of models and observations of He I lines came from Porter 
et al. [TT] who compared model predictions with measured fluxes for 100 lines with n < 20 from 
the Orion nebula. They found an average difference of 6.5% between all observed and predicted 
lines and 3.8% difference for the 22 most accurately measured lines. 

Zhang et al. [12\ [T3] present a method to diagnose temperatures from He I line ratios. 
They write parametric forms for He I line ratios vs. temperature and conclude that the best 
diagnostic is obtained from the A7281/A6678 ratio. The results of the diagnostics have an 
intrinsic uncertainty of ~1000 K at around 10 000 K owing to optical depth effects on the 
He I A3889 (2s 3 S -> 3p 3 P°) line. When comparing the results from the A7281/A6678 ratio 
with those from A7281/A5876 a systematic error of (~500 K) arises, which may be related to 
the uncertainty on the optical depth of the A3889 line, or perhaps from some other source. 
Interestingly, in a sample of 48 PNe the helium temperature is systematically lower than that 
obtained from the hydrogen H/? decrement by an average of 4000 K. 

4. Second and third row elements (C-Ar) 

Density and temperature diagnostics in PNe are commonly based on line ratios among 
collisionally excited lines. The best known density diagnostics are [S n] A6717/A6731 and [O n] 
A3729/A3726 that are sensitive to densities of the order of 10 3 cm~ 3 , [CI III] A5517/A5537 for 
N e ~ 10 4 cm~ 3 , and [Ar iv] A4711/A4740 for N e ~ 10 4 - 10 5 cm- 3 . 

The [O n] ratio that arises from transitions among the 4 £3/2 ground level and the 2 D 5 / 2 
and 2 D 3 / 2 levels was a subject of controversy for the last few years. It was commonly assumed 
that L5*-coupling was a valid approximation for terms of the ground configuration of the + 
system; thus the ratio of collision strengths from the ground levels to the 2 D 5 / 2 an d 2 D 3 / 2 levels 
was given by the statistical weights, i.e. 1.5. However, McLaughlin & Bell [14] claimed that 
relativistic effects enhanced the ratio of collision strengths to 1.93, with profound effects on N e 
determinations for low surface brightness H n regions. Such claims were contested by evidence 
from the extensive literature survey of Copetti &: Writzl [15J and the observational campaign 
of Wang et al. [16]. The issue seems settled now with a new calculation by Pradhan et al. 
|17] , that accounts for all dominant relativistic effects and confirms the earlier predictions of the 
LS-coupling approximation. The spectroscopic survey by Wang et al. also served to test the 
A- values for dipole forbidden transitions of [O n]. They found that the transition probabilities 
of Zeippen [18] best fit the observations, the differences being of only a few percent, while the 
results of later calculations by the same author ([E]) and by W. Wenaker [20J look problematic. 
This latter set is what is currently available through the NIST database (|21j). 

Wang et al. also compared the electron densities obtained from various line ratios and 
obtained very good agreement between N e ([0 n]), iV e ([S n]), and N e ([Cl in] ) . By contrast 
A e ([Ar iv]) values yield densities systematically higher than those from the other diagnostics, 
which sheds some doubts on the accuracy of the iV e ([Ar iv]) A- values. 



5. The iron-peak elements 

Fe and iron-peak elements are important constituents of gaseous nebulae. Depending on the 
excitation of the nebula, iron is frequently seen in stages from Fe I to Fe vn. The observed ions 
of nickel span about the same range. Other less abundant elements of the group are occasionally 
identified as well. An extreme example of rich spectra is the r/ Carinae ejecta (see article by 
Gull in this volume). 

On this subject the contributions of Sveneric Johansson could hardly be overstated. He has 
been the driver of extensive and detailed research of energy levels and transition rates for the low 
ionization stages of iron-peak species. Data often are taken for granted as they becomes easily 
available through various databases, yet their determination from laboratory work requires very 
skillful people and substantial funding. Johansson's publications on the subject are too numerous 
to list here, but some of the most important contributions are: on the spectra of Sc n [22J, Fe I 
[23], Fe II |24j . Ti II [25], Co II [26], and on the determination of /-values, lifetimes, and radiative 
rates for forbidden transitions (e.g. for Fe II [23 [25] US ED], Ti II [31], Ni I [35]). Johansson has 
also done seminal work the understanding of the Lya excitation mechanism of Fe n and other 
iron-peak species in astronomical spectra (e.g. [33 ], 134" ] 135 ] l5o]). 

Spectra of singly ionized iron-peak species are particularly complex as they respond to a 
variety of excitation mechanisms. These are H Lya fluorescence, continuum fluorescence, self- 
fluorescence by overlapping Fe n transitions, and collisional transitions among low and high 
levels. In moderately dense plasmas (iV e ~ 10 7 cm~ 3 ) the populations of high levels involved in 
fluorescence are redistributed by collisionally induced transitions through highly excited pseudo- 
metastable levels [37]. Furthermore, proper modeling of Fe n spectra requires large systems, 
approaching one thousand levels and complete and accurate atomic data. This has led to a 
diverse set of models with little consensus among them and limited understanding about their 
accuracy. For ions of atoms other than Fe and Ni the first detailed spectral models are now 
being created, e.g. [38j [39]. A good part of this work has been done under the auspices of the 
IRON Project [10]. A summary of the data is presented in Fig. 1. 

More work is also in progress by various groups like the IRON Project, the FERRUM Project 
at Lund University, at NIST, and the group at Queen's University Belfast. 

Another modeling challenge in dealing with dense spectra of iron-peak species is the treatment 
of the transfer of line radiation from the point of emission to the boundary of the atmosphere. For 
iron-peak species, in addition to continuum opacity and resonant scattering, one must consider 
fluorescence by line overlaps, such as the well known H Lya effect in Fe n [33]. While writing 
the present review, I have estimated the number of overlaps between absorption lines arising 
from low excitation levels, typically populated under nebular conditions, of Fe ii-iv and lines 
from H I, He II, O n and O III, dominant species in typical nebular spectra. Line wavelengths 
awere taken from the NIST database and do not represent the entire spectra. For line widths 
of 20 km/s the number of overlaps are 592, 243, and 49 for Fe II, III, and iv respectively, while 
for widths or 60 km/s the overlaps are 1722, 762, 524 for the same ions. This is the range of 
line widths in nebular spectra. The huge number of line overlaps demonstrates the need for a 
detailed treatment of radiative transfer of Fe spectra. The case also applies to other iron-peak 
species. 

In modeling spectra one also needs to solve for ionization equilibrium of each element. For that 
reason, detailed level specific photoionization cross sections, that account for complex resonance 
structures, are being computed for both ground and excited states using the R-matrix method. 
Cross sections are available for Fe I-IV [50] EH [52] [53] and Ni II [53] . Recently, much progress has 
been made at the Instituto Venezolano de Investigaciones Cientificas (IVIC) on the calculation 
of cross sections for neutral Sc, Ti, and Cr. 




Figure 1. Spectral models for low ionization iron-peak species 



6. Hyperfine-induced spectra 

Observations of hyperfine-induced transitions in PNe have been of great interest in recent years. 
This is because these spectra allow us to determine relative isotopic compositions to be compared 
with predictions of nucleosynthesis. This, in turn, enables us to disentangle enrichment by stellar 
nucleosynthesis from the primordial composition in observed present total abundances. For 
example, the abundance of 3 He in our Galaxy could be used to test big bang nucleosynthesis 
and constrain the baryonic density of the Universe. However, the evolution of 3 He in the 
Galaxy is a longstanding open problem due to the observed discrepancy (by two orders of 
magnitude) between theoretical yields of low-mass stars and the measured abundances in H n 
regions [551 156| . A proposed solution of this problem could be in suppressing, by non-standard 
mixing mechanisms, the production of 3 He during the red giant branch and/or asymptotic giant 
branch phases of stars of mass up to ~2 Mq \57\ 158] . An observable consequence of such an 
scenario is that the ratio 12 C/ 13 C in the ejecta of PNe should be much lower than in the standard 
case. For a 1 Mq star, the predicted ratio is ~5, in contrast with the standard ratio of 25-30 
[59]. 

Isotopic abundances can be derived from beryllium-like ions through observations of the 
hyperfine-induced transition 2s2p 3 Pq — 2s 2 1 Sq within the UV0.01 multiplet. This transition 
can only occur through coupling of the electronic total momentum with nuclei of nonzero spin. 
This happens in the spectra of 13 C (with nuclear spin 1/2) but not in 12 C. Thus, 12 C/ 13 C 
abundance ratios can be derived from spectra showing the hyperfine transitions together with 
the stronger 2s2p 3 P 2 ° - 2s 2 (M2) and 2s2p 3 P 1 ° - 2s 2 1 S (IC) transitions. Palla et al. 
[60] observed these lines in NGC 3242 using STIS on board the Hubble Space Telescope. From 
this, they derived 12 C/ 13 C =38. The same lines from the isoelectronic ion N iv were measured 
by Brage et al. [61] from a STIS spectrum of NGC 3918. In this case both stable isotopes 



of nitrogen ( 14 N and 15 N) have non-zero nuclear spin. Although no abundances were derived, 
this work served to confirm the theoretical A-value for the hyperfine transitions. More recently, 
Rubin et al. [62] analyzed spectra of PNe in the IUE archive and determined 12 C/ 13 C = 4.4±1.2 
for one object and established lower limits to this isotope ratio in 23 objects. 

Another study of hyperfine induced transitions, but this time in Ne-like Al iv, was reported 
by Casassus et al. |63| . The two stable isotopes, 27 A1 and 26 Al, have nuclear spins 5/2 and 5 
respectively. Thus, hyperfine structure causes the 3 P2 - 3 Pi transition at 3.66/iin to split into 9 
components, 5 of which could be resolved spectroscopically. From the observations, Casassus et 
al. calculated an upper limit to the 26 A1/ 27 A1 abundance ratio of 1/33 in NGC 6302, which is 
in good agreement with the predicted value of 1/37 from stellar evolution models. 

7. Neutron capture elements 

One of the frontiers of spectroscopy of PNe is the study of n-capture elements presumably 
synthesized by the AGB progenitor of the nebula. The progress on this area since the seminal 
work of Pequignot &; Baluteau [65] has been slow owing to the need for high spectral resolution 
(> 10 000) and the absence of atomic data and models. Despite the difficulties, Dinerstein |65] 
was able to identify lines of [Kr III] and [Se iv] in the K-band of the near infrared spectra of 
NGC 7027 and IC 5117. Soon after, Dinerstein & Geballe [66J also detected [Zn IV] lines in the 
near infrared spectrum of NGC 7027 . To date these lines have been identified in almost one 
hundred PNe [571 [68]. 

On the theoretical side, the calculation of collisional data for heavy species represents an 
important challenge because intermediate coupling representations are inappropriate in many 
cases, at the same time as relativistic effects become too large to treat in the Breit-Pauli 
representation. Thus, researchers must turn to fully relativistic codes in the Dirac formalism. 
Recently, Badnell et al. [HU] put the Dirac Coulomb R-matrix package of Berrington et al. [70] 
on the same footing as the traditional Breit-Pauli R-matrix codes developed by the IRON/RmaX 
team. Nonetheless, these sort of calculations are still far from routine work, owing to the size 
of the atomic representations in JK -coupling. 

8. Additional considerations 

This review would not be complete without pointing out some important general issues of 
relevance to the present subject as well as recently developed tools for atomic spectroscopy. 

The accuracy of recombination rate coefficients was in recent years an important source of 
concern for modelers for their effects on calculations of ionic fractions. The issue, however, 
seems to have reached a point of stability as different theoretical methods, e.g. the unified R- 
matrix approach of Nahar and Pradhan and the Thomas- Fermi-Dirac approach in the version of 
Badnell and collaborators [71] . begin to yield consistent results and to release large amounts of 
data. Such theoretical results also seem to agree with recent experimental determinations (e.g. 
|72j). A warning must be raised though regarding low temperature dielectronic recombination, 
that is dominated by near threshold autoionizing channels, whose positions are difficult to get 
accurately from theory. More experimental work on this area is needed to benchmark current 
and forthcoming data (e.g. [73]). 

On its origins, one of the explicit objectives of the IRON Project was to provide a complete and 
reliable set of data for positive ions of interest in nebular infrared spectroscopy. A comprehensive 
review of these data is presented by Badnell et al. [75]. The entire dataset produced by the 
IRON Project will soon be available through TIPbase |75| . 

A very useful tool for nebular spectroscopy, the EMILI package for emission line identification, 
has been developed by [75]. This tool should expedite the analysis of high spectral resolution 
and signal-to-noise spectra, and enable the identification of faint poorly known features. 



9. Conclusions 

The question of what fractions of Lya photons are self-absorbed and removed by dust in 
photoionized regions is rather important, as the treatment of hydrogen spectra in photoionization 
models is one of the main limitations to their accuracy. There are also various questions in regard 
to heating and photo-evaporation of grains in PNe. Further research on this is fundamental for 
entering a new level of detailed understanding of PNe and H n regions. 

Temperature diagnostics with H and He are very important as they yield new information 
on the long lasting problem of temperature fluctuations in gaseous nebulae. It is important, 
though, to sort out the apparent inconsistencies between temperature diagnostics from different 
recombination series of lines or line ratios. Moreover, it seems plausible that temperature 
fluctuations in nebulae could be accompanied by density variations. Thus, it is important 
that both quantities be derived simultaneously from the same diagnostics. 

As atomic models have become available for essentially all of the most prominent species in 
spectra of gaseous nebulae it is important to check that diagnostics of temperature and density 
out of various species all yield consistent results. This would allow us to determine whether 
the atomic data have reached ultimate accuracy. A flag of warning must be raised, however, 
in regard to neutral species and ions with ionization potentials below that of hydrogen because 
their spectra could be affected by photo-excitation by continuum radiation. Such is the case of 
N I [77] that Copetti & Writzl [15] tried to analyze on the same grounds as ionized species. 

Iron and iron-peak ions are as important as ever, particularly Fe II, which is prominent 
throughout astronomy, beyond PNe research. These ions, however still pose great challenges. 
The open 3d structure of these species is difficult to describe owing to various factors: (1) large 
numbers of strongly correlated LS terms shortly spaced in energy from the ground state; (2) 
strong angular correlations among configurations such as 4s 2 , 4p 2 , 4d 2 , 4p4f that lead to very 
demanding computations; (3) strong core-valence interactions in the vicinity of excitations of 
3p electrons onto the unoccupied 3d orbital, leading to the so-called giant 3p— >3d resonance; 
(4) electronic correlations among 3d electrons such as the radial distribution of electrons in a 
3d^ +1 configuration differs from that of electron in a Sd^nl. More experimental work is needed, 
particularly in the measurement of cross sections that guide theoretical work (e.g. |86j). 

The importance to modern astronomy of studying isotopic abundances and n-capture 
elements can hardly be overstressed. Yet, such work imposes extreme demands on spectroscopic 
instruments and astronomers. Carefully coordinated efforts between atomic physicists, 
astronomers, and instrument designers would be desirable on this subject. 
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